Ejection dynamics of a ring polymer out of a nanochannel 



Junfang Shcng^ and Kaifu Luo^'B 

CAS Key Laboratory of Soft Matter Chemistry, 
Department of Polymer Science and Engineering, 
University of Science and Technology of China, Hefei, Anhui Province 230026, P. R. China 

(Dated: September 20, 2011) 

We investigate the ejection dynamics of a ring polymer out of a cylindrical nanochannel using both 
theoretical analysis and three dimensional Langevin dynamics simulations. The ejection dynamics 
for ring polymers shows two regimes like for linear polymers, depending on the relative length of the 
chain compared with the channel. For long chains with length A'^ larger than the critical chain length 
Nc, at which the chain just fully occupies the nanochannel, the ejection for ring polymers is faster 
compared with linear chains of identical length due to a larger entropic pulling force; while for short 
chains [N < Nc), it takes longer time for ring polymers to eject out of the channel due to a longer 
distance to be diffused to reach the exit of the channel before experiencing the entropic pulling 
force. These results can help understand many biological processes, such as bacterial chromosome 
segregation. 

PACS numbers: 87.15.-v, 82.35.Lr, 87.15.H- 



I. INTRODUCTION 

The properties of a polymer confined in a nanochannel 
have attracted broad interest [11-01 because they are of 
fundamental relevance in polymer physics and are also 
related to many biological processes, such as double- 
stranded DNA genomes packaging inside the phage cap- 
sid @ , polymers transport through nanoporejp, [l3| and 
viruses injecting their DNA into a host cell \Vi\ . 

The importance of cyclic structures in biological 
macromolccular science is strikingly demonstrated by the 
existence of circular DNA, cyclic peptides and cyclic 
oligosaccharides and polysaccharides [IJl. Ring closure 
of a polymer is one of the important factors influenc- 
ing its statistical mechanical properties. Understand- 
ing the static and dynamic properties of ring polymer 
is a challenging problem due to the difficulties inherent 
to a systematic theoretical analysis of such objects con- 
strained to a unique topology. The scaling behavior of 
isolated, highly diluted, ring polymers has been studied, 
des Cloizeaux [l^l, Deutsch 3] and Grosberg [l^ dis- 
cussed the effect of topological constraints on the proper- 
ties of ring polymers, and found that the topological con- 
straint and the excluded volume have similar effects. The 
radius of gyration for large single ring polymers obey the 
same scaling relationship as that of linear chains [ij, [l^ , 
although this is not true for ring polymers in a melt or 
ring polymer brushes [l6l - [l8j . 

Ring closure acts as an important role in a wide range 
of biophysical contexts where DNA is constrained: segre- 
gation of the compacted circular genome of some bacteria 
18 1 , formation of chromosomal territories in cell nuclei 
20|, compaction and ejection of the knotted DNA of a 
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virus [2l|, |22|, migration of a circular DNA in an elec- 
trophoresis gel |23| or in a nanochannel [24 1. 

After three decades of intensive research, the confor- 
mational properties of a self-avoiding polymer chain con- 
fined in a slit or in a c ylin drical nanochannel are relatively 
well understood. |25l429| . However, a deeper understand- 
ing of the basic properties of ring poly mer in confined en- 
vironments is a field in its infancy [30, llll ■ Only few stud- 
ies have addressed semiflexible ring polymers. Ostermeir 
et al. [s^l investigated the internal structure of semi- 
flexible ring polymers in weak spherical confinement and 
found buckling and a conformational transition to a fig- 
ure eight form. Fritsche and Hecrmann [ssj examined the 
conformational properties of a semiflexible ring polymer 
confined to different geometrical constraints and found 
that the geometry of confinement plays a important role 
in shaping the spatial organization of polymers. Most 
recently, we have found the helix chain conformation of 
flexible ring polymers conflned to a cylindrical nanochan- 
nel, and demonstrated that the longitudinal size along 
the channel for a ring polymer scales as Na{a/D)^^^, 
the same as that for a linear chain but with different 
prcfactors. Here D is the radius of the channel, N the 
chain length and a the Kuhn length of the chain [s^ . 
We further gives the theoretical ratio value 0.561 of the 
longitudinal size for a ring polymer and a linear chain of 
the same TV. 

As to the dynamics of the polymer under confinements, 
Milchev et al. [131 have investigated the ejection of linear 
chain out of nanopore using Monte carlo simulation and 
found that the ejection dynamics depends on the chain 
length. Unlike its linear polymer counterpart, the dy- 
namics of confined ring polymers is still lacking, although 
many bimolecules are circular. To this end, in this work 
we study the ejection dynamics of a ring polymer con- 
fined in a nanochannel by means of analytical techniques 
and Langevin dynamics simulations. The basic questions 
associated with this process are the following: (a) what's 
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FIG. 1: A schematic representation of the initial and the finial 
conformation of a ring polymer ejection from a cylindrical 
channel in three dimensions. Here, the channel diameter D = 
5, the channel height h = 20.5 and the polymer length A*' = 
100. 



effect of the chain length and the channel length on the 
ejection dynamics? (b) what's the difference of the ejec- 
tion dynamics for ring polymers compared with the linear 
one? For a fixed channel, which one is faster compared a 
ring polymer with a linear chain of the identical length? 

We believe that this work is interesting and important 
for understanding biological systems with more complex- 
ity, such as viruses injecting their DNA into a host cell, 
the behavior of DNA inside phages or the spatial organi- 
zation of the bacterial nucleoid in E. coli. 



II. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard- Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction between 
beads is modeled by a short range repulsive LJ potential: 
UL.j{r) = 4£[(f - i^ f] + e for r < and for 

r > 2i/V. Here, a is the diameter of a bead, and e 
is the depth of the potential. The connectivity between 
neighboring beads is modeled as a FENE spring with 
UpENsir) = — ^fci?Qln(l — r'^/Rl), where r is the dis- 
tance between consecutive beads, k is the spring constant 
and Rq is the maximum allowed separation between con- 
nected beads. 

We consider a schematic representation as shown in 
Fig. [U where a ring polymer is confined in a cylindri- 
cal channel with one end sealed. The nanochannel and 
the sealed surface are described by stationary particles 
within distance a from one another which interact with 
the beads by the repulsive Lennard- Jones potential. The 



particle positions of the nanochannel and the sealed sur- 
face are not changed in the simulations. 

In the Langevin dynamics simulation, each bead is sub- 
jected to conservative, frictional, and random forces, re- 
spectively, with [H mfj = -y{ULj+UFENE)-^'^i+Ff'- 
Here m is the bead's mass, ^ is the friction coefficient, 
Vi is the bead's velocity, and is the random force 
which satisfies the fluctuation-dissipation theorem. In 
the present work, the LJ parameters e, a, and m fix the 
system energy, length and mass units respectively, lead- 
ing to the corresponding time scale Ilj = (mcr^/e)^/^ 
and force scale e/tr, which are of the order of ps and pN, 
respectively. The dimensionlcss parameters in the model 
are then chosen to be Rq ~ 1.5, fc = 15, ^ = 0.7. 

In our model, each bead corresponds to a Kuhn length 
(twice of the persistence length) of a polymer. For a 
single-stranded DNA (ssDNA), the persistence length of 
the ssDNA is sequence and solvent dependent and varies 
in a wide range, to our knowledge, usually from about 
1 to 4 nm. We assume the value of a ~ 2.8 nm for a 
ssDNA containing approximately four nucleotide bases. 
The average mass of a base in DNA is about 312 amu, so 
the bead mass m w 1248 amu. We set ksT ~ 1.2e, which 
means that the interaction strength e is 3.39 x 10"^"'^ J 
at actual temperature 295 K. This leads to a time scale 
of 69.2 ps and a force scale of 1.2 pN. The Langevin 
equation is then integrated in time by a method described 
by Ermak and Buckholz (36j . 

We initially fix the last monomer of the linear chain but 
anyone of the ring polymer at the sealed bottom of the 
nanochannel, while the remaining monomers are under 
thermal collisions described by the Langevin thermostat 
to obtain an equilibrium configuration. In order to learn 
the mechanism of chain ejection out of the nanochannel, 
the link of the monomer with the bottom of the chan- 
nel is removed, then the chain is released to diffuse along 
the channel. The residence time r is measured, once all 
monomers pass the opening at x = h and leave the chan- 
nel. Typically, we average our data over 700 independent 
runs. 



III. RESULTS AND DISCUSSION 

A. Scaling arguments 

1. Longitudinal size of a polymer in infinitely long 
nanochannel 

According to the blob picture, for a linear polymer con- 
fined in a infinitely long three-dimensional nanochannel 
of diameter D, the chain will extend along the channel 
axis forming a string of blobs of size D. The center of 
the blob is on the axis of the nanochannel. For each 
blob, D = Ag'^a due to the dominant excluded volume 
effects, where g is the number of monomers in a blob, a 
is the Kuhn length of the chain, u is the Flory exponent 
in three dimensions, and A is a constant. Thus, each 
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blob contains g = i'^)" monomers, and the number of 
blobs is ni, = N/g = N{^)i. The free energy cost 
for the chain confinement is proportional to the num- 
ber of blobs, thus the free energy in units of ksT is 
J" = BiN{Aa/Dy/'', with Bi being a constant. The blob 
picture then predicts the longitudinal size of the linear 
chain to be = mD = ND{^)i = {Aa)iND^-i. 
Using = 3/5 in three dimensions, we obtain 



Rni = {Aa)^ND- 



(1) 



In order to model the the chain conformation for a ring 
polymer confined in a nanochannel, we have extended the 
blob picture [3^ . For a ring polymer, the chain will ex- 
tend along the channel axis forming two strings of blobs 
of _D/2,the two strings of blobs show helix structure. For 
each blob of size D/2, D/2 — Ag'^a with gr being the 
the number of monomers in a blob. Here, the same pref- 
actor A for ring polymers and linear chains is due to 
the same solution environment. Thus, each blob con- 
tains gr = (2^)" monomers, and the number of blobs 
is nil = N/g.,, = N{^^)i. The free energy cost in units 
of fcsT is = BrN{^^)i , with Br being a constant. 

By geometrical analysis, the distance between two suc- 
cessive layers is ^D, and so the total length occupied by 
blobs in the channel is R- — y^rti n 

V2 W I 



When — is very large 
1^ — 3/5 in three dimensions, we obtain 
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Therefore, the longitudinal size along the channel for 
a ring polymer scales as _R|| ~ Na{a / D)'^/'^ , the same as 
that for a linear chain but with a different prefactor. The 
ratio of the longitudinal sizes along the nanochannel (or 
the prefactors) for a ring polymer and a linear chain is 



R\\ 



^^2t =0.561. 



(3) 



If using more accurate value of = 0.588, we have = 

-"-II, i 

0.575. 

The simulation results [U confirm the above predic- 
tions and give {AuY^'^ = 1.367±0.009 for the parameters 
used in the model. 



Ejection dynamics of a polymer confined in a 
nanochannel 



Intuitively, for the ejection of a polymer out of a 
nanochannel, the dynamics is controlled by the relative 
length of polymer compared with the channel height h. 



There exists a critical polymer length Nc-, where the poly- 
mer just fully occupies the channel, namely {Nc) = h. 
Thus, for linear chains the critical length Nc^i is 



Nr. 
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{Aa)iD 

while for ring polymers the critical length Nc^r reads 
4\/2/i 1.782/1 



Nrr = — 



2's{Aa)iD-i {AaY^D- 



(4) 



(5) 



The ratio of the critical length for the ring polymers and 
the linear polymer is 



Nrr 4^2 



^ = 1.782. 
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(6) 



Short chains with chain length oi N < Nc are initially 
fully confined in the nanochannel while long polymers 
with chain length oi N > Nc initially occupy the whole 
channel with several segments outside the channel exit. 
For long chains with N > Nc, the ejection is a driven 
process where the pulling force / is from the entropy 
and is induced by already ejected monomers @. For 
short chains with N < Nc, polymer needs to move to the 
channel exit by a diffusive process, and then experiences 
a pulling force as for long chains. 

We assume the ejection process to be quasi- 
equilibrium. For long chains, the pulling force can be 
estimated from the free energy of a chain partially con- 
fined in the nanochannel with the innermost monomer 
being at distance x from the channel exit. For long lin- 
ear chains, x = n{t){Aay''" D^^^^'^ and the free energy 
J" = Bin{t){Aa/Dy^''kBT = Bi^x, with n{t) being 
the number of monomers inside the channel at time t. 
The differential of the free energy allows an estimate of 
the pulling force 



fi=B^ 



D 



(7) 



It is worthy of noting that / is independent of the tail 
length as well as h, but inversely proportional to D. 
For long ring chains, x = n{t) {Aa)i D^~i and the 
free energy J" = Brn{t){2Aa/DY/''kBT = iy^Br^x. 
Thus, the pulling force 



fr =4V2Br- 



D 
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We further have the ratio of the pulling force for long 
ring polymers and linear chains 



f_r_ ^ AV2Br 
fl 



Bi 



(9) 



which is only determined by the universal prefactors for 
ring polymers and linear chains. 



During the ejection process, the pulling force induced 
by the tail is balanced by the total friction. Namely, for 
long chains we have 



(10) 



where ^ is the friction coefficient per monomer. Taking 
into account the relationship of x{t) and we obtain 
the ejection time 



for long linear chains, and 



Tlo 



(11) 
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for long ring polymers. Therefore, the ratio of the ejec- 
tion time for long ring polymers and linear chains is 



'^long.r 



^ = 0.315—, 13) 

'long, I r>r^ -^r 



where = 3/5 is used. 

As noted above, for short polymers {N < Nc), it un- 
dergoes a diffusive process before the first segment exiting 
the channel, and subsequently the ejection process driven 
by a pulling force. Accordingly, we divide the total ejec- 
tion time T into two parts, Ti for the diffusive process 
and T2 for the driven process. 

For the the diffusive process of short linear chains, ti 

is 

jh-Rli)^ m[h-{A<j)iND'-i]^ 
'^1,1 = — ;rr; = ttt^ ' (14) 



2kBT 



with Ddiff = ^^^j- being the diffusion constant. In addi- 
tion, for the driven process can be written as 



T2 J 



2Bi{AGfl-kBT 



^N^{Aa)iD^-i 
2BikBT 



(15) 



Here, T2^i is negligible compared to ri,; for quite short 
chains, and then the ejection time t; « ti^/. Based 
on the differential of the ejection time with N, ^ = 

[h - (Acr)^iVL>i-^][/i - 3{Aa)iND^-i]=0, we obtain 
the critical chain length Nc,i = as shown in 

Eq. and another resolution 



(16) 



3{A<j)iyD^-i 

where the residence time r/ reaches to its maximum value 

2Ch^ 



27{Aa)iD^'ikBT' 



(17) 



For the the diffusive process of short ring chains, ri is 

_ (/i-i?ll,,)2 _ N^[h- ^2i{A<j)iyND^-i]^ 
^^'^ ~ 2Dd,ff ~ 2fc^T ' 

In addition, for the driven process T2 can be written as 



(18) 
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T2,r = 



_ ^N^{2Act)i:D 
2Br{2Aay/''kBT ^ QABrkBT 



(19) 



Again, T2,r is negligible compared to ti_,. for quite 
short chains, and then Tr ~ ri^^- Based on the dif- 
ferential of the ejection time with iV, ~ [h — 
^{2Aa)i ND^-i][h - ^{2Aa)i ND^-i] = 0, we ob- 
tain the critical chain length Nc,r 
Eq. ([5]) and another resolution 



as in 



4:V2h 



(20) 



3{2Aa)iD^~y 
where the ejection time t reaches to its maximum value 



27{2Aa)iD^-ikBT 



Thus, we have 



n: 



4^2 



^ 1.782. 



(21) 



(22) 



B. Simulation results 

The average ejection time r as a function of the ring 
polymer length N for different channel diameters {D — 
5, 7 and 9) at fixed channel height h = 20.5 and for differ- 
ent channel heights {h = 20.5,30.5, and 40.5) at channel 
diameter D = 7 are shown in Fig. [2^ and Fig. [2]3, respec- 
tively. The two pictures show that ejection time increases 
with the increase of channel diameter and channel height. 
Moreover, we get a special polymer length N* at which 
the ejection time meets its maximum. Fig. [3] shows the 
plot of TV* against hD'^/^ for different D and h. All the 
data points collapse on the same line, which is in agree- 
ment with the prediction in Eq. ([20)) . The line plotted 
in Fig. [4] proves the prediction in Eq. (|2T|) . 

As noted before, there exists a critical polymer length 
Nc at which the polymer just fully occupies the channel. 
Short chains (A^ < Nc) are initially fully confined in the 
nanochannel while long polymers [N > Nc) initially oc- 
cupy the whole channel with several segments outside the 
channel exit. From the platforms in Fig. [21 we obtain 
the ejection time Tiong for long polymers {N > Nc). 

Fig. [5^ and Fig. [SJd show the scaling plot of Tiong with 
ring polymers and linear chains, respec- 
tively. For different polymer lengths, channel heights 
and channel diameters, all data points collapse on the 
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FIG. 2: (a) Average ejection time r as a function of the poly- 
mer length for the channel height h = 20.5 and different 
channel diameters, (b) The same for channel diameter D = 7 
and different channel heights. The full lines are plotted using 
Eqs. (fT8)) and (fT9)) . 




100 125 



FIG. 3: Scaling of Nc at which the ejection time owns its 
maximum as a function of hD^^^ for different h and D. Here, 
h = 20.5, 30.5, and 40.5 and D = 5, 7 and 9. 
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FIG. 4: Scaling of the maximum of ejection time of different 
polymer sizes as a function of h?' D^^^ for different h and D. 
Here, we use D = 5, 7 and 9 for h = 30.5 and 40.5, and D — 4, 
5, 7 and 9 for /i = 20.5. 
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FIG. 5: The ejection time nong as a function of h?D^^^: (a) 
ring chain with iV = 250, 300 (D = 5, 7 and 9 for /i = 20.5, 
30.5 and 40.5); (b) linear chain with N = 250, 300 (D = 3, 5 
and 7 for ft = 20.5 and 30.5), N = 350, and 360 (D = 9 for 
ft = 40.5). 
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FIG. 6: Average ejection time as a function of the polymer 
length for the linear chain and the ring polymer. Here, D — 7, 
and h = 20.5. The full line plotted to fit the data points for 
linear chain is based on Eqs. (fT4)) and (fTS)) . 
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FIG. 8: Squared distance x'^{t) of the last monomer from the 
channel exit during the ejection as a function of the time for 
the ring polymer and the linear chain and different channel 
diameters. Here, h = 20.5 and A'^ = 250. 
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FIG. 7: The force / exerted on the innermost monomer of the 
ring polymer and the linear chain. Here, A'^ = 300, D = 3, 
and h = 20.5. 



same line in Fig. [S^ and Fig. [SId, respectively. These 
results confirm the predictions in Eqs. (|TT|) and (|T2|) . 
In addition, the slopes are 0.042 and 0.053 for ring 
polymer and linear chain, respectively. This indicates 
= M|2 ^ 0.792. Based on Eqs. (HH) and ([12]), we 

have Br = 1.60, Bi = 4.02 and thus ^ = 0.398 using 
the parameters C = 0.7, T = 1.2 and {Aa)^^" = 1.367. 
Moreover, we further obtain ^ = 2.250 through Eq. ©, 
which demonstrates that the driving force induced by 
confinement for long ring polymers is larger than that 
for linear chains. Using Eqs. ((T5)) and to fit curves 
in Fig. [21 we find that the numerical results are qualita- 
tively described by theoretical findings. 

To compare the ejection dynamics for ring polymers 
with that for linear chains, we show the ejection time as 
a function of the chain length N for D = 7 and h = 



20.5 in Fig. [B] One does see characteristic differences: 
for short chains {N < Nc), it takes longer time for ring 
polymers to eject out of the channel than that for linear 
chains; while for long chains {N > Nc), linear chains 
need longer time. These findings are in agreement with 
the predictions in Eqs. ([T^ and (P^ . Ring polymers has 
smaller than that for linear chains of the same N and 
thus ring polymers must diffuse longer distance to reach 
the exit of the channel. When the chain length is larger 
than the critical chain length {N > Nc), the force exerted 
on the residual segments for ring polymer is larger than 
that for linear chain due to the smaller blob size in the 
channel for ring polymers than that for linear chains as 
predicted ^ = 2.250. The platform of the force at small 
time t shown in Fig. [7] for both ring polymer and linear 
chain confirms this prediction. 

In Fig. we find ^ = II = 1.417, = 

= 1.456, which are predicted as 1.782. The dif- 
ference may be from the non-equilibrium process of the 
ejection. In addition, ~ ^ ~ 1.719 and it is pre- 
dicted as 1.782 in Eq. ® . 

The mean-squared distance a;^ {t) of the last monomer 
with respect to the channel exit against elapsed time af- 
ter the release of the last monomer is shown in Fig. [5] 
for both ring polymers and linear chains. The lines plot- 
ted according to the curves are based on the equation 

x^{t) = a:^(0) — ^ ^ 1 for a linear polymer 

and the equation x'^{t) = a;^(0) — /■■(^^^'^^^^^ f foj- 

a ring polymer, which indicates that the Eq. pop can 
correctly describe the ejection dynamics. 

In order to know the details in the ejection process, 
we record the number of residual monomers inside the 
channel in the total process, n{t)/n{0) (normalized by 
its value at t ~ 0). We see that the short ring polymer 
experiences a diffusion process before it starts to eject 
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FIG. 9: Numbers of monomers n{t) inside the channel nor- 
malized by its value at t = as a function of the time for 
(a) ring polymers with different lengths; (b) the ring polymer 
and the linear chain of length A*' = 100. Here, D = 3 and 
h = 20.5. 



out of the channel, corresponding to the platform in the 
plot for = 30 as shown in Fig. [5^. When the chain 
length N > Nc, the ejection process is faster for the ring 
polymer than that for a linear one, which can also be 
inferred from the portion of residual monomers at time 
t, as presented in Fig. 

Fig. [TOk shows the histograms of the ejection time for 
ring polymers with different chain lengths. The ejection 
time distribution for polymer of length = 30 has a 
long tail and is much wider than that for = 50. The 
ejection time distributions for ring polymers and linear 
chains at both short and long chain lengths are given in 
Fig. [TUb and Fig. [TUb . respectively. For short chain 
N = 20, it takes longer time for ring polymer to leave 
the channel than that for linear chain, and the ejection 
time distribution for the ring polymer is wider and has a 
long tail. For long chains N = 300, however, the result 
is opposite, reflecting the larger driving force for the ring 
polymer than that for the linear chain. 

Nature not only imposes geometrical constraints on 
biopolymers by confinement through cell membrane, the 
cell nucleus or viral capsid, but also exploits the advan- 
tages of certain underlying chain topologies, such as the 
ring structure. In fact, E. coli has a rod-shaped geometry 
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FIG. 10: The distribution of the ejection time: (a) ring poly- 
mers of different chain lengths; (b) the ring polymer and the 
linear chain of length A*' = 20; (c) the ring polymer and the 
linear chain of length A'' = 300. Here, D = ?, and h — 20.5. 



and its chromosome is not a linear polymer but a circular 
one. Based on Monte Carlo simulations, Jun and Mul- 
der [l^ addressed a basic physical issue associated with 
bacterial chromosome segregation in rod-shaped cell-like 
geometry. By simulations of two ring polymers, in the 
same setting as the linear ones and they found that two 
ring polymers segregate more readily than linear ones in 
confinement. According to our above theoretical analy- 
sis and simulation results, for ring polymers confined in 
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a cylindrical nanochannel the blob size for ring polymers 
is smaller than that for linear polymers, which indicates 
that during the chromosome segregation the driving force 
for ring polymers is larger than that for linear one, lead- 
ing to faster segregation. 

IV. CONCLUSIONS 

We investigate the ejection dynamics of a ring polymer 
out of a cylindrical nanochannel using both theoretical 
analysis and three dimensional Langevin dynamics sim- 
ulations. The ejection dynamics for ring polymers shows 
two regimes like for linear polymers, depending on the 
relative length of the chain compared with the channel. 
For long chains with length larger than the critical 
chain length Nc, at which the chain just fully occupies 
the nanochannel, the ejection for ring polymers is faster 
compared with linear chains of identical length due to a 
larger entropic pulling force; while for short (N < Nc), 
it takes longer time for ring polymers to eject out of the 
channel due to a longer distance to be diffused to reach 



the exit of the channel before experiencing the entropic 
pulling force. These results can help understand many 
biological processes. 

Our results should enable a new understanding of the 
conformational statistics and dynamics of confined ring 
biopolymers such as DNA. The concrete graph about 
ring polymer confined in a nanochannel needs more deep 
study so as to realize many complex problems in both 
biochemistry and theoretical study. Our findings are of 
interest for (molecular) biological/biochemical, technol- 
ogy as well as physics reasons. 
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